Real space tomography of the primordial Universe with cluster polarization 
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We describe how a survey of the polarization of the cosmic microwave background induced by 
Compton scattering in galaxy clusters can be used to make a full spatial reconstruction of the 
primordial [z ~ 1089) matter distribution inside our surface of last scattering. This "polarization 
tomography" can yield a spatial map of the initial state of the Universe just as gravitational collapse 
was beginning to drive structure formation. We present a transparent method and simple formulas 
from which one can compute the 3D primordial map in real and in Fourier space, given a 3D map 
of the polarization due to galaxy clusters. The advantage of the real space reconstruction is that it 
is free from the statistical uncertainties which are inherent in the Fourier space reconstruction. We 
discuss how noise, partial sky covering and depth of the survey can affect the results. 

PACS numbers: 98.65.Cw, 98.70.Vc, 98.80.-k, 98.80.Es 



Introduction The latest observations of the cosmic mi- 
crowave background (CMB) P, SJ and of the galaxy dis- 
tribution in the low-redshift Universe [H, Q| have shown a 
spectacular agreement with the standard theory of struc- 
ture formation @|. In this scenario, the initial condi- 
tions were set up by inflation, which created a nearly 
scale-invariant spectrum of Gaussian density perturba- 
tions with tiny variance (a ~ 10~ 5 ), and all the visible 
structure we see today developed through gravitational 
instability from the inflationary seeds. 

However, these two sets of observations are made at 
two very different moments in time, and they probe very 
different regions of space: while the CMB is formed es- 
sentially through the Sachs- Wolfe effect at the time of 
recombination (z ~ 1089), galaxies on the other hand are 
observed only in relatively recent times (z < 5.) More- 
over, the CMB gives us a picture of the thin spherical 
shell known as the last scattering surface (LSS), which 
corresponds to the original location, at the time of recom- 
bination, of the CMB photons we observe on Earth today, 
while the galaxy distribution, in practice, can only be ob- 
served well inside our LSS. This fact leads to a statistical 
"leap of faith" in the argument for the standard scenario: 
we believe that the matter distribution inside our LSS is 
similar to the matter distribution on our LSS, and there- 
fore a typical configuration for the matter distribution at 
z ~ 1089 would lead to a map of the large-scale structure 
at z < 5 which has similar statistical properties to the 
one we actually observe. This leap can be a particularly 
long one if the CMB anisotropies, which provide a picture 
of the matter distribution on our LSS, shows suspicious 
anomalies @, 1 S ES El El El, 0, El, HE El El • 

Kamionkowski and Loeb (2(| were the first to real- 
ize that it is possible, in principle, to recover informa- 
tion about the primordial matter distribution inside our 
LSS using cluster polarization of the CMB - see also 
[II 11 SH, 11, SI[. This is because the CMB photons, as 



'Electronic address: abramo@fma.if.usp.br 



they Compton-scatter on free electrons in a galaxy clus- 
ter (the Sunyaev-Zel'dovich effect - see [1ES3|), acquire 
a polarization which reflects the quadrupole of the tem- 
perature distribution of the incident photons (28| . Hence, 
observing the polarization pattern in a galaxy cluster is 
tantamount to observing the quadrupole of the CMB in 
the LSS seen by that cluster. In principle, this means 
that by sampling all the LSS's of clusters inside our LSS 
through cmb polarization, we could reconstruct the tem- 
perature distribution of photons in the whole volume in- 
side our LSS, and thus reconstruct the matter distribu- 
tion in that volume - see Fig. 1. 

Notice that a similar argument could be made for the 
CMB polarization which is induced at the time of last 
scattering through the local temperature quadrupole as 
well as for the polarization from the epoch of reionization. 
However, the former is highly degenerate with the infor- 
mation we already have from the CMB temperature itself 
(see [Hj] for a reconstruction scheme based on the CMB 
temperature), and the latter is a large-scale effect which 
is widely spread over the line of sight (as the ionization 
depth in this case is only significant over horizon-scale 
distances) , making it difficult to trace the polarized pho- 
tons back to a particular spacetime event over the past 
light-cone. Nevertheless, if reionization is patchy and 
inhomogeneous, polarization combined with 21-cm ob- 
servations can be used much in the same way as cluster 
polarization to reconstruct the primordial matter distri- 
bution [3CJ | - and formulas very similar to the ones we 
show here can be used in that case as well. 

In practice, the reconstruction scheme based on clus- 
ter polarization has been difficult to implement for sev- 
eral reasons. First and foremost, because the amplitude 
of the polarized signal is at most a few percent of the 
temperature fluctuation signal for typical cluster opti- 
cal depths [H, SI SI, SI . Second, foregrounds such as 
lensing (34[, extended galactic emissions (3f|SES3, clus- 
ter peculiar velocities [38l . SI and filamentary structures 
[3ll ] . can be hard to extricate from the data. And third, 
the methods for recovering the matter distribution from 
polarization maps lacked transparency [H SH, S3 ■ 
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We believe that these difficulties will be overcome. 
First, the latest technological advances have made it 
possible to start searching for the B-mode of polar- 
ization induced by gravitational waves from inflation 
[40L l4ll l42l [43| - which will be more than enough sensitiv- 
ity to detect polarization at the levels found in clusters. 
Optic and X-ray data, which are crucial to determine the 
optical depths and redshifts of the clusters, should make 
de-lensing significantly easier at the arcminute scales we 
are interested 0, |45j • Filamentary structures can be 
avoided by cutting out regions of the sky with very high 
gas columns or projected surface density of clusters. Sim- 
ilarly, galactic foregrounds can be either directly sub- 
tracted with multi-frequency observations, or a mask can 
be applied onto the maps. Cluster peculiar velocities can 
also be dealt with in multi-frequency surveys [22^ . 

As for the reconstruction of the matter distribution, we 
will present a transparent method and simple formulas 
for the reconstruction of the matter distribution both in 
Fourier as well as in real space. We have considered the 
contributions from the Sachs- Wolfe and the integrated 
Sachs- Wolfe effects (SWe and ISWe, respectively). It is 
trivial to consider surveys with partial sky covering, as 
well as any type of redshift binning one may employ. We 
also study the effect of adding noise to the polarization 
data and its impact on the reconstructed maps. 

Cluster polarization Since only the Stokes parameters 
Q and U are relevant in cosmology [H, [47], [Hf , we define 
the dimensionless complex polarization of a cluster at 
redshift z in the direction n as: 



P(z,h) 



Q(z, h) —iU(z, h) 

T(z)T z ,n 



(1) 



where T{z) = 2.726(1 + z) K and t z ^ is the optical depth 
for the cluster at redshift z and position h. 

At linear order in the cosmological perturbations (we 
neglect the Riess-Sciama effect), the polarization for a 
given cluster depends on the gravitational potential $ 
on the cluster's LSS through the SWe, and it depends on 
the fluctuations of the potential inside the cluster's LSS 
through the ISWe - see Fig. 1. 

Usually one describes the gravitational potential in 
terms of its Fourier spectrum: 



$[r(z, n)} — 



d 3 k 
(2tt) 3 / 2 



-ik-r 



(2) 



The Fourier transform of the gravitational potential has 
Gaussian random P.D.F.'s set up by inflation, with 
($(fc)$*(fc')) = k~ 3 P(k)S(k-k'), and P{k) = Ak 11 ^ 1 a 
nearly scale-invariant function. 

For this problem it is very useful to describe the grav- 
itational potential in the volume inside our LSS (which 
includes the LSS's of all clusters) in the following way: 



$(r)=$^m(r)Y* ro (n) , 



(3) 
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FIG. 1: Last scattering surfaces seen by an observer at the ori- 
gin at n = 770 (solid black circle) and by two clusters at times 
771 and 772 (dashed grey circles - dashed red and blue circles 
in the online version.) The variable x is physical distance in 
units of the distance to our LSS, x — (77 — r/ rec )/(rio — r] rec ). 
The distance from the origin to the LSS of a given cluster at 
x spans the interval |1 — 2x\ < x' < 1. 



where r — rjo — n{z) and 77 is conformal time. The coeffi- 
cients b£ m (r) are well-behaved functions of the radius if 
is a smooth function. Since the radius of our LSS, 
at z rcc ~ 1089, is i?LSS = ^0 — *7recj it is useful to write 
all physical distances in terms of this maximal radius, 
x = r/i?Lss so x lies in the interval < x < 1. 

We can write &(k) in the same way as Eq. ©, in terms 
of coefficients bg m (k). The formal relationship between 
btm{x) and bp m {k) is a Hankel transform: 



bt m {k) = \j— (— i) / dxx je(kx) b lm (x) . (4) 



The inverse relation is found by exchanging k «-> x 
and multiplying by (— 1) . The fact that $(r) is real 
translates into the conditions b* lm = (—l) e+m be _ m and 

A decomposition analogous to Eq. ([3]) can be made for 
the cluster polarization as well, except that here the spin- 
2 character of the polarization field demands that we use 
the spin- weighted spherical harmonics instead [251 ] : 



(5) 



In terms of the coefficients pi m (x) and be m (k), the clus- 
ter polarization takes on a very simple expression: 



pe m (x) = / dkKe(x,k)be m (k) , 



(6) 



where the Fourier space Kernel Ki(x, k) is given by: 



K t (x,k) =i f 60Trd?-k 2 f e {kx)A 2 (x 1 k) , (7) 



£.771 
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with [25|: 
fi(kx) 



~ E (-l) (A ^ )/2 (2A + l) (8) 



x=e-2,£,e+2 
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and 



A 2 (z,fc) = J2[*(l-*)] - 6 



1 j ii ■ r; c « m dG(x") 
dx j 2 [k{x - x)\ ^ 



(9) 

The first and second terms in Eq. © stand for, re- 
spectively, the SWe and the ISWe. The growth function 
G(x) = D[z(x)]/a[z(x)} is normalized to unity today and 
its derivative is non-negligible only in relatively recent 
times (z <j 5) as dark energy becomes more relevant. 

In fact, it turns out that the sum in Eq. ([8]) is a particu- 
lar realization of the spherical Bcsscl function's recursion 
relations, and the result is: 



fi(kx) 



(£ + 2)1 n (kx) 



6(£-2)\ 3Q(kx) 2 



(10) 



This window function, which mediates the exchange of 
power between the "orbit" angular momenta £ > 2 and 
the spin angular momentum of the polarization field, is 
identical to the one found in Ref. [24[ for the correlation 
function of the polarization coefficients, and it ensures 
that the power is conserved, Y^T=2^ + l)/ 2 = 1- It is 
also the same (up to a normalization factor) as the radial 

function ef 1 of [48| , which shows that the cluster-induced 
polarization is a pure E (or gradient) mode. 

Polarization tomography in real space Even per- 
fect knowledge of the polarization field on the whole vol- 
ume of our LSS does not eliminate the statistical errors 
in the Fourier space coefficients b( m (k) - as Eq. (01 is 
cut-off at x = 1. These uncertainties are a 3D analog of 
the 2D "cosmic variance" that plagues CMB anisotropies. 
However, even though these statistical uncertainties can 
be large, especially for the longest-wavelength modes, we 
can still reconstruct the 3D spatial map of the primordial 
fluctuations inside the LSS with as much accuracy as the 
observations allow. This is a nontrivial statement, as the 
polarization we observe from a cluster carries only the in- 
formation of the temperature quadrupole as seen by that 
cluster, projected along the line of sight. Nevertheless, 
it turns out that this is just enough information to allow 
(at least in theory) for a complete spatial reconstruction. 

Plugging Eq. ([4]) into Eq. ([6]) and using a generaliza- 
tion of the Weber-Shafheitlin integral [49[ we obtain: 

Plm(x) = / dx' K e (x,x')be m (x') , (11) 
Jo 

where the real space kernel Kg(x, x') is given by an SWe 
piece and an ISWe piece. It turns out that Ki is real, 



which is another manifestation of the fact that cluster 
polarization is a pure E mode. 
The SWe kernel is: 



K, 



(SW) 



(£ + 2)1 



,.' 3 



2-k(£-2)\ x(1-x) 3 

x 6(\x' + x\ - |1 - a:|)] 8(\1 -x\- \x' - x\) 

x sin 2 i/;P £ ~ 2 (cosV') , (12) 

where P e n is the generalized Legendre polynomial, and 
cosV> = [x 2 + x 12 — (1 — x) 2 ]/ (2xx r ) is the cosine of the 
angle between x and x' in a triangle with sides x, x' and 
1 — x - see Fig. 1. The step- functions in Eq. (fT2")l . which 
vanish unless 1 — 2x\ < x' < 1, automatically imply 
causality and ensure that the polarization in each point 
can only depend, through the SWe, on the gravitational 
potential on the LSS of that point. 
The ISWe part is given by: 



K 



(ISW) 



-Sir 2 



dx" - 



1 



(£ + 2)1 x' 
2tx(£-2)\ ~x J x (x"-xf 

x 9(\x' + x\- \x" - x\)] 9(\x" -x\- \x' - x\) 
dG(x") 



dx 



sin 2 V'i^ (cosV') , 



(13) 



where cos-0' = [x 2 +x l2 ~ (x" — x) 2 ]/(2xx') is the cosine of 
the angle between x and x' in the triangle with sides x, x' 
and x" — x. As was the case for the SWe kernel above, it 
is straightforward to show from the step-functions in Eq. 
(fT3"|) that the ISWe kernel vanishes if x' > 1. Detailed 
geometric interpretations of the real space kernels can be 
found in a forthcoming publication [501 ] - 

Notice that by virtue of the step-functions in the ker- 



nels above, the variables . 



and x" all lie in the interval 



[0, 1], so both p£ m and bt m have compact support. This 
means that the linear problem of finding the coefficients 
b£ m (x') given the data pe m (x) is completely well-defined 
in real space. This is in contrast with the problem of 
finding the Fourier space coefficients bi m given the same 
data: because the p^ TO 's in Eq. ^ can only be inte- 
grated over the interval < x < 1, the inversion of Eq. 
© is singular. Evidently, this is a manifestation of the 
statistical uncertainties of the 6^ m 's. 

Eqs. (fTTj) - (fT3)) are the main results of this paper. Next 
we discuss applications of these formulas. 

Map reconstruction In practice, the integrals ([6]) and 
(|lip become sums over bins, so we have to solve a linear 
problem of the type P tm (i) = £\ K t (ij)B tm y) , with K t 
some matrix that need not be square if we want to com- 
pute less parameters Bg m than there are data points 
Pim(i)- But even if the matrix Ki is square, in general 
it is singular. This is a severe problem for the Fourier 
space kernel, which has many singular eigenvalues due to 
the statistical uncertainties in the long-wavelength modes 
A -Rlss L 50j . For the real space kernel we find that 
there are essentially no singular eigenvalues, which al- 
lows, at least in principle, for an exact inversion. 
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t = 10 | 




a = 0.50 



0.1 Vo 
z 

FIG. 2: Noisy estimator Kf{± + N t )Ki for £ = 10 (top) 
and £ = 30 (bottom), assuming zero errors (black stars) and 
assuming Gaussian errors with variances ao = 0.05, 0.2, 0.5 
(boxes of dark, medium and light grey.) We used 500 realiza- 
tions of Nt for each ao, and we use / s k y = 2/3. In this example 
we have assumed a hypothetical survey with 100 redshift bins 
up to 2 = 10, and we recover information on 50 redshift bins. 

For linear problems which involve rectangular or sin- 
gular square matrices, we can employ the pseudo-inverse 
[51|. The pseudo -inverse K + of a matrix K reduces to 
the inverse in the case of square non-singular matrices, 
and it gives the least-squares solution B = K + ■ P to 
the linear problem P = K ■ B. In this framework it is 
trivial to consider noisy data: the observations become 
Pim — > Pf m = (t + Ne)-Pe m , where Nt is the noise matrix 
which is diagonal if the errors are uncorrelated. 

In the presence of noisy data the estimator Bg m be- 
comes a noisy estimator Bf m , which can be related to 
the actual parameters Bg m that we want to measure by: 

B? m = K+ (1 + Ni)Ki ■ B em . (14) 

Here Ni is a Gaussian random variable with zero mean 
and variance given by: 

WsWb( w -* + ^ V2 ' (15) 

where ao includes factors such as instrument noise, num- 
ber of pixels and redshift, and / s k y < 1 is the fraction 
of the sky mapped by that instrument. The term / s k y 
inside the parenthesis in Eq. (fl~5|) subtracts cosmic vari- 
ance, which is not relevant for the accuracy of the actual 
map reconstruction. 

Since for the real space kernel K + K is nearly identical 
to the unit matrix (unlike the case of the Fourier space 



kernel) , one can easily find that the covariance matrix of 
the estimator B^ m is given by: 

Cov i B hn)i0 — < J 1$ij'^ B t7nkB* lmk . (16) 

k 

One can recognize the sum over redshift bins in the r.h.s. 
of Eq. (fl6|) as the projected Cg - i.e., J dzCe(z). 

If the matrix Kg~(t + Ng)Kg is close to a unit matrix, 
then the noisy estimator will be a good approximation 
to the true parameters. We can evaluate the goodness 
of the noisy estimator by looking at the diagonal of that 
matrix. 

We have done this for all the modes between I = 2 and 
I = 50, using errors ao = 0.05, 0.2, 0.5, in a hypothetical 
survey up to z = 10 with 100 bins, and in which we re- 
construct the bim's in 50 redshift bins. In this example, 
for 1 = 1 the reconstruction fails, mainly because of poor 
sky coverage. However, as i increases sky coverage be- 
comes less of a problem and the reconstruction becomes 
progressively better. Fig. 2 shows the accuracy of the re- 
construction for I = 10 and I — 30. A deeper survey does 
not significantly improve the reconstruction, but a survey 
shallower than z 3 impacts negatively the reconstruc- 
tion at low redshifts, as the LSS's of low-rcdshift clusters 
are entirely at high redshifts. The severest limitation to 
the reconstruction scheme is sky coverage. 

Conclusions We have shown how to reconstruct the 
primordial density field inside our LSS using cluster po- 
larization. In principle, a complete reconstruction is pos- 
sible with this method. In practice, sky coverage and 
noise limit the accuracy of the reconstruction, especially 
for low I's and for intermediate redshifts (z ~ 1). Nev- 
ertheless, a reconstruction of the primordial density field 
in the local Universe (z < 0.5) and on angular scales 
2° ^5 < 20°, is feasible with the next generation of po- 
larimeters which are being used to search for the B mode 
from gravitational waves. 

The primordial field is highly correlated with the 
present matter distribution. Direct comparison of the 
two is a key test of our models of structure formation. 
In particular, correlating the primordial map with the 
present locations of clusters and superclusters tests the 
evolution of structures, as potential wells only move due 
to nonlinear effects. Finally, cross-correlating the primor- 
dial map with tracers of the growth of structure is a new 
test that constrains the parameters of dark energy. 
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